home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48hor2 / gausselm.doc < prev    next >
Text File  |  1995-03-31  |  1KB  |  32 lines
  1. GAUSSELM - Gaussian elimination and QR factorization 
  2. Author: Robert Brunner (brunner@uirvld.csl.uiuc.edu) 
  3. Date: 7 Oct 91 02:39:23 GMT 
  4.  
  5. The following directory contains two useful linear algebra routines and three 
  6. subroutine functions. 
  7.   
  8. QR:  QR decomposes a matrix on the stack, A, into an orthogonal matrix Q and a 
  9. upper triangular matrix R such that A=QR, using Gram-Schmidt orthogonalization. 
  10. The program only works for full-rank matricies. 
  11.   
  12. LU:  LU decomposes the matrix from the stack, A, into a permutation matrix P, a 
  13. lower triangular matrix L, and an upper triangular matrix U, such that PA=LU. 
  14. Pivoting is only performed when necessary, and the first non-zero pivot is 
  15. selected, rather than always choosing the largest pivot.  By doing this, 
  16. pivoting is avoided when not needed, which is useful for solving small problems 
  17. given for homework in linear algebra classes.  The method does exhibit poor 
  18. numerical stability, so the program is not practical for large matrices. 
  19.   
  20. DCMP:  DCMP decomposes a matrix [[row 1][row 2]...[row n]] into a list of 
  21. vectors {[row 1][row 2]...[row n]}.  Doing this speeds up the other two 
  22. routines. 
  23.   
  24. RCMP:  RCMP reverses DCMP. 
  25.   
  26. SWRW:  SWRW switches the two specified rows of a matrix decomposed by DCMP. 
  27.   
  28. 'GAUSSELM' BYTES gives checksum #33398d, 1458 bytes 
  29.  
  30. Drop me a line if you find this useful. 
  31. brunner@uirvld.csl.uiuc.edu 
  32.